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Abstract 

A larger than expected forward-backward asymmetry in rapidity is observed in top quark pairs 
produced in proton-antiproton collisions at the Tevatron. The asymmetry is seen in both the top 
quark distribution A t FB and in the distribution of charged leptons A FB from top quark decay. 
In this paper, we study the kinematic and dynamic aspects of the tight relationship of the two 
observables arising from the spin correlation between the charged lepton and the top quark with 
different polarization states. We also consider two benchmark new physics models, an axigluon 
model and a flavor-changing W' model. These models could explain the values of both A FB and 
A FB . We emphasize the value of both measurements, and we conclude that a model which produces 
more right-handed than left-handed top quarks is favored by the present data. 
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I. INTRODUCTION 



The observed forward-backward asymmetry in rapidity A FB of top quarks jl, 2] at the 
Fermilab Tevatron deviates by about two standard deviations (2a) from standard model 
(SM) expectations [3]. After corrections for detector acceptance and resolution, A FB in 
the center-of-mass (cm.) frame is 15.8 ± 7.5% at CDF and is 19.6 ± 6.5% at DO Q. 
In addition to A FB , the DO group also reports a positive forward-backward asymmetry of 
charged leptons from top quark decays of A FB = (15.2 ± 4.0)% compared with the small 
value 2.1 ±0.1% from simulations of the SM [2]. The definition of A FB is 



.4 



Tip — 71b 



FB 



(1) 



71 F + n B 

where n F {n B ) is the number of events with qiyi > (qiyi < 0), and qi and are the 
sign and rapidity respectively of the charged lepton from the semileptonic decay of a top or 
anti-top quark in the ti production. 

In this paper, we investigate the kinematic and dynamic relationship between the two 
observables A FB and A FB . We study quantitatively the influence of the top-quark boost 
on the kinematics of the charged lepton, showing how the distribution of leptons in the 
laboratory frame is related to the polarization state of the top quark parent. We show 
that current data on the ratio of the two asymmetries favor models in which more right- 
handed than left-handed top quarks are produced. The fact that A FB , A FB , and the ratio 
A FB /A FB are larger than the SM predictions indicates that the charged lepton strongly 
prefers to move in the same direction as the top quark from which it originates. This result 
can arise if right-handed top quarks % |5j play a significant role in A FB or if a non-standard 
mechanism produces more highly boosted top quarks at the Tevatron, as we explain below. 

Many new physics (NP) models have been proposed to explain the enhancement of A FB , 



such as flavor- changing Z' [6j, W \L\ and axigluon G' |8MlO| models. 1 The first two mod- 
els produce predominantly right-handed top quarks, whereas the axigluon model generates 
unpolarized top-quarks. It is important to validate these models at the Large Hadron Col- 
lider (LHC) and/or at the Tevatron. For example, the heavy flavor-changing Z' mt) 



1 The next-to-leading order quantum chromodynamics corrections to the process of qq — > tt induced by the 
flavor-changing Z' and W are calculated in Ref. ll| and Ref. 12(, respectively, with the result that the 



NP prediction at the leading order is reliable. 
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avored because it predicts too much same-sign top quark pair production at the 



14|. In this paper, we focus on how consistently the NP models can describe both 



^fb and A FB . 

We begin in Sec. [Ill with a discussion of the angular distribution of decay leptons, first 
in the rest frame of the top quark and then after the top quark is boosted in rapidity and 
transverse momentum. We pay particular attention to left /right polarization state of the top 
quark because the final distribution of leptons in the laboratory frame, after the top quark 
is boosted, depends significantly on the top quark's polarization state. In Sec. IIIH we derive 
the relationship of the lepton asymmetry A FB and the top quark asymmetry A FB separately 
for the left- and right-handed polarization states of the top quark. Different models of new 
physics produce top quarks with different proportions of left- and right-handed polarization. 
We use two such models, an axigluon model and a W model, in Sec. [TV] to deduce their 
different expectations for the ratio of the lepton and top quark asymmetries. Our conclusions 
appear in Sec. |Vl We emphasize the value of making measurements of both A FB and A FB 
because their ratio can be related through top quark polarization to the underlying dynamics 
of top quark production. 

II. KINEMATICS 

The charged lepton in top quark decay is a powerful analyzer of the polarization of the 



top quark 15j . In the rest frame of a top quark, the distribution in the polar angle #hei of a 



decay lepton £ + is 

1 dT _ 1 + X t cos fl he i ^ 

T d COS #hel 2 

where A t denotes the top quark helicity. Here, X t = + is for a right-handed top quark (t#) 
while Xt = — for a left-handed top quark (t^). The angle is measured with resect to the 
direction of motion of the top quark in the overall center-of-mass system of the ti production 
process. The distributions are shown in Fig. []Ja). The charged lepton from a right-handed 
top quark decay prefers to move along the top quark direction of motion, while a lepton 
from a left-handed top quark moves preferentially against the top quark direction of motion. 
In the rest frame of the top quark, 75% (25%) of charged leptons from tn (t£) decay follow 
the top quark direction of motion, i.e. cos^hei > 0. 

Once the top quark is boosted, the angular distribution of the charged lepton relative to 
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FIG. 1: (a) cos#hei distribution in the top quark rest frame for both tz and £r. (b) cos 9a 
distribution in the boosted frame for a top quark with E t = 200 GeV. 

the direction of motion of the top quark is sensitive to the energy of the top quark E t (or 
equivalently its velocity (3). We derive 

dT l-/3cos0t« + At(cos0tf-/9) 



Td cos d tl 2 7 2 (1 -/3 cos 9 



(3) 



where (3 = yl — mf/Ef, 7 = E t /m t . In Fig. [TJb) we plot the distribution in cos 6^ of 
the charged lepton, where the angle cos9 t e is the angle between the charged lepton and its 
parent top quark in the boosted frame. As an illustration, fixing the energy of the top quark 
to Et = 200 GeV, we find that about 60% of £ + follow the top quark for a ti, and almost 
100% for a t R . 

The distribution of charged leptons in the laboratory frame depends on the top quark 
kinematics, including the top quark energy and its rapidity, and the top quark polarization. 
The probability for finding a positive charged lepton in the forward region when it originates 
from a top quark with a velocity (3, rapidity yt, and polarization At is defined as 

Aft 

where N l F (Ng) denotes the number of leptons i in the forward (backward) region in the 
laboratory. After lengthy algebra, it can be shown that the ratio R F is 
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where 
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(6) 



To illustrate the effect of the top quark boost, we plot in Fig. [2] the fraction Rp as a 
function of y t . We choose two characteristic top quark energies, E t = 200 GeV and 600 GeV. 
The former energy represents top quarks produced around the threshold region, while the 
latter pertains for highly boosted top quarks. Note that ?/ max = 0.53 for E t = 200 GeV. 
When a top quark moves perpendicular to the beam line, i.e. y t = 0, there is an equal 
number of leptons in the forward and backward regions, leading to Rf = 0.5, independent 
of E t and the polarization of the top quark. 

For right-handed top quarks tn, Rf increases rapidly with y t in the region of y t > 
because most of the leptons move close to the direction of motion of the top quark after 
being boosted to the lab frame; this result is shown by the black solid lines in Fig. [TJ We can 
also see that when E t becomes larger, i.e. the top quark is more energetic and the lepton is 
more boosted, Rf rapidly reaches its maximum value 1. 
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FIG. 3: The ratio of the charged lepton in the forward and backward region as a function of the 
top quark rapidity for top quarks with fixed transverse momentum pr = 10, 50, 100, 300 GeV. For 
a fixed pt = 50 GeV, the figures show that around the region of yt ~ 0.2, the fraction of charged 
leptons in the forward region is about 75% for a right-handed top quark while 45% for a left-handed 
top quark. 

On the contrary, in the case of t^'s, the ratio Rf does not vary significantly with y t owing 
to the anti-boost effect on the charged lepton. For E t = 200 GeV, the boost causes charged 
leptons to distribute nearly uniformly, and Rp is around 0.5, as seen in the red-dotted curve 
in Fig. [2](a). When the energy of t^'s is large enough, the large boost forces most of the 
charged leptons from top quark decays to move along the top quark direction of motion, 
even if they move against the top quark direction of motion in the top quark rest frame. 
The boost yields a large value Rp in the region of large y t , as shown by the red-dotted curve 
in Fig. Et^b). The competing influences leave the curve slightly below the t^j curve. 

In Fig. |3l we show how R F varies with p l T and y t . The distributions for right-handed 
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top quarks £r's do not vary greatly with p T because most of the charged leptons follow £r. 
However, the shapes of the curves for left-handed top quarks, which are the focus in the 
discussion below, are very different between the low p T and high p T regions, as is seen in the 
red-dotted lines. As the top quark moves forward, i.e. yt > for fixed p T , the boost becomes 
more significant as the energy of the top quark is increased. Therefore, more leptons are 
forced to move along the direction of the top quark. On the other hand, some fraction of 
the decay leptons which are initially in the forward/backward region (ye > 0/ye < 0) will 
then be in the backward/forward region by definition. In summary, there are two factors 
which affect Rf'- the boost and the rearrangement of the distribution of charged leptons 
in the forward (ye > 0) and backward (ye < 0) region. The former always increases Rf 
while the latter may increase or decrease the Rf depending on how energetic the top quark 
is at y t = 0. Generally speaking, when the boost is not significant (low p T and small y t ), 
Rf decreases when y t increases from y t = 0, as we can see in the drop in the red-dotted 
curves in Fig. [3|a) and[3|b). When the boost is big enough, Rf always increases with y t . 
The platform-like behavior around yt — in Fig. [3](c) arises because the leptons accumulate 
nearly uniformly around the axis of motion of the top quark when p T = m t /\/3 ~ 100 GeV. 
Therefore the ratio Rp is rather stable as the top quark changes its direction of motion 
direction around y t = 0. 



III. A FB AND A FB 

The observed positive top-quark asymmetry A FB indicates more top quarks are produced 
in the forward region than in the backward region of rapidity. Both tn and ti can generate 
a positive lepton asymmetry A FB . However, as shown in Fig. [2], tp would need a large boost 
along the proton beam line (i.e. in the large forward rapidity region) to overcome the fact 
that most of the charged leptons from its decay move against it in its rest frame. A right- 
handed top quark t# can yield a positive A FB even for top quarks near the ti threshold region. 
Therefore, the large positive top quark and lepton asymmetries A FB and A l FB observed by 
the DO collaboration indicate that the top quark polarization may be playing a non-trivial 
role. In this section we present a general analysis of the correlation between A FB and A e FB , 
to prepare for a better understanding of the numerical results derived from NP models to 
be shown in Sec. IIVI 
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The top quark asymmetry A FB can be expressed as a sum of contributions from the SM 
and NP as: 

#FB = ^FB P >< R + X (1 - J2), (7) 

where 

sm _ iVl M ~ iVl M ,t, np _ iVg - NT p _ <t P ^ 

" jySM + jySM FB " ^NP + ^NP > « jySM + jyNP > W 

with -/VJyL and NpT™ being the numbers of events in which the top quark moves with 
y t > 0(y t < 0) in the SM and induced by NP, respectively, and N^^ NF ^ is the total number of 
events predicted in the SM (induced by NP). The NLO QCD contribution to the production 
process qq — > ti could generate a value A FB ~ 5%, which is much less than the central 
value of experimental data. 

To somewhat simplify the discussion of the correlation between A t FB and A e FB , we assume 
in this section that A FB is generated completely by NP, but all SM contributions (including 
the NLO QCD effects) are retained in the numerical calculations we present. 

The contributions to A FB from different polarizations of top quarks can be separated as: 



where 



rt-FB 



At, NP 



A t L , NP . A t R , NP] w d / n \ 

Pt L A fb + Pt R A£ B x R, (9) 



A-FB 



' PA t = ^np • (10) 

NP iv tot 



V£ + 

Here, A FB KP denotes the forward-backward asymmetry of the top quark with polarization X t 
generated only by NP, while p\ t is the fraction of top quarks with polarization At in tt events 
induced by NP. One advantage of decomposing A FB into different top quark polarizations is 
to monitor the chirality of the couplings of NP particles to top quarks. Another advantage 
is to make the connection between A FB and A FB more transparent. 

As discussed in Sec. [II] the ratio R F depends on the top quark kinematics (/3, y t and A t ). 
To compute the probability for a charged lepton in the forward region, one must convolute 
the top quark production cross section with R F on an event-by-event basis, i.e. 

Ar« g, tf* = j N*(J3, y t , \ t )R F Xt ((3, yt) d$, (11) 

where N tt labels the ti production rate for a top quark with specific kinematics (/3, y t , Xt) 
and $ stands for the phase space. The lepton asymmetry A l FB generated by a top quark 
with polarization X t is 
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Np ® Rp Xt + N% g> i^ At - iVp ® i^ At - N Xt (g> i^ At 



NP 



A^< + JV** 



NP 



A^ ® (W/ 4 - l) + N% ® (2i$ At - 1 



(N Xt -N Xt )® (2Rp Xt - 1 



NP 



N Xt + N Xt 



Here, 



NP 



l-i4 At (/3, m ), 



and we use the following relation between Rp Xt and Rg Xt in our derivation, 



R% Xt (y t )=R F ^(-y t ). 



(12) 



(13) 



(14) 



The quantities iVj^ and ATg in the convolutions in Eq. (TT2]) should be understood as 
the distributions N u ({3,y t , A t )0 (y t ) and N tt ((3,y t , \ t )Q (—y t ), respectively, where (x) 
is the Heaviside step function. The quantity N F — Ng should be understood as 
[N*(j3,y t , \t) - N a (f3, -y t , A*)] 6 (y t ). Because R e F Xt in Eq. flT2]) cannot exceed 1, we have 
< ^fb- When R F Xt is close to a constant Tic, e.g. T^c ~ 1/2 around the tt threshold 
(E t ~ 200GeV) for left-handed top quark or 1Z C ~ 1 for a highly boosted top quark, Eq. (fT2"|) 
can be simplified as 
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(15) 



Equation (fl5|) and Fig. [2] show that: 

• A% ~ when the tt pair is produced around the threshold region; 

• Ap 1 ^ < Ap§ ~ v4^ B in the large region. 

Although Eq. (fT5j) is approximate, it helps in understanding the NP prediction obtained 
from a complete numerical calculation. 
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IV. NEW PHYSICS MODELS: AXIGLUON AND W 



In this section we focus on two models of new physics, an axigluon model |8MlO| and 
a flavor- changing W model [7J . We examine how these NP models can accommodate the 
values of both A FB and A FB measured by the DO collaboration. 

In the axigluon (G') model we assume for simplicity that the interaction of the axigluon 
to the SM quarks is purely pseudo-vector-like and can be written as 

C = g s {& qi^q + g h QYisQ) g;, (16) 

where q denotes the first two generation quarks in the SM and Q the third generation quarks. 
The coupling g s is the usual strong coupling strength; gi and gh are the coupling strength 
(normalized to the QCD strong coupling g s ) of the axigluon to the light quark (q) and the 
heavy quark (Q), respectively. 

The helicity amplitudes of the processes qq — >• g — >• tt and qq — > G' — > tt are written as 
M g (X q , Xg, X t , Af), and Ma'(X q , Xq, X t , Af), where A, = + represents the right-handed helicity 
of particle i and Aj = — the left-handed helicity. The total helicity amplitude is 

-gWh) 



M(X q ,Xq, At, A*) =gXaP. 



2+A +A 

cd 



•Mg(X q , Xq, X t , At) + 2 — r~ Fi — Mg'{\, Aq, Af, At 

S — m QI + ITflQlY Qt 



,(17) 



where tfj is the generator of the color SU (3) group; a, b, c and d are the color indexes of 
q, q, t and t, respectively. The non-vanishing helicity amplitudes are 



M g (- + —) = -Mg(+-++) = y/l-f3 2 sm6, 

Mg(+ ) = -Mg(- + + + ) = Vl-P 2 Sin 9, 

Mg(~ + - + ) = Mg(+ ~ +-) = "(I + COS#), 

Mg(-++~) = Mg(+--+) = (1-COS0), (18) 



and 



M G ,(+ - +-) = M G >(- + -+) = (3(1 + cos9), 

M G >(- + +-) = M G i+ - -+) = 0(1 - cos9), (19) 



where = \/l — Ami I $ an d ^ i s the polar angle of the top quark in the cm. frame of the 
tt pair measured relative to the initial state quark. 
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FIG. 4: Correlation between A FB and A FB for (a) the axigluon and (b) the W' models. The point 
corresponding to the DO data is also shown. The numbers within the parentheses label the lower 
and upper limits of the mass of the NP object. The statistics for the ratio of predicted A FB to 
A FB for the G' and W models are shown in (c). For comparison, the SM values are A l FB ~ 5% 
(off the left side of the plots in (a) and (b), and A FB ~ 2%. 

The absence of deviation from the SM expectation in the measured m# distribution 
indicates the axigluon should be heavy and broad. The axigluon's contribution to tt 
production is therefore through interference with the SM channel. The interference effect 
becomes largest in the region of large m t t, i.e. ~ 1. Therefore, the last two equations 
of Eq. (fT8|) dominate. When < ffJc, the denominator of the axigluon propagator is 
negative, and the square of the interference term in the overall amplitude is proportional to 

s 



[2 9l g h {l + cos#) 2 - 2g l9h (l - cosfl) 2 ] 



m 



(20) 



G> 



The term linear in cos 6* is 4gig h s cos 9/ (s — m 2 G ,). The product g\g h must be negative to 
obtain a positive A FB js-10]. 

The top quarks are generated unpolarized owing to the pseudo-vector coupling of the 
axigluon to the SM fermions, and 



Pt L 



PtR 



At 

A t L , NP _ At R , NP _ A FB 
^FB ~ ^FB ~ £, 



FB 



> 0. 



(21) 



Since the tt cross section is greatest near the threshold region where A F l B ~ and A F l B 
A FB , the expression for A l FB becomes 



4b ~ pt L 



A% NP (2TZ C -l)xR + PtR A% NP (2TZ C - 1) x R 



-A 



FB- 



(22) 
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We plot our axigluon model predictions for A FB and A FB in Fig. IU(a). We first scan the 
theoretical parameter space (gi, gh and mc) to fit Tevatron data on A FB and the ti total 
production cross section within 1 a. These parameters are then used to calculate A e FB . The 
figure shows a clear correlation between A FB and A e FB . The best fit to the correlation is 

A FB ~ 0.47 x A FB + 0.25% . (23) 

To fit both A FB and A l FB within la, the mass of the G' must be greater than 1 TeV. For 
masses this great, top quarks from G' decays are highly boosted and cause most of the 
charged leptons to move along the direction of the top quarks. We remark here that if the 
G' is found as a resonance in the ti mass distribution, the chirality structure of its coupling 
to ti can possibly be determined at the LHC [if]]. 

A different class of NP models to explain the ti forward-backward asymmetry is based on 
t-channel kinematics. Such models involve large flavor- changing interactions. A model with 
a non-universal massive neutral vector boson Z' 6] is one of the possibilities. However, it is 
disfavored because it implies a large rate for same-sign top quark production at the 7 TeV 



LHC 



13], not supported by data [14]. 



We consider in this paper a flavor- changing W which couples an incident <i-quark to the 
produced t-quark 

C = gigRd-fPntWl + h.c. , (24) 

where g<i is the weak coupling. In the W model, in addition to the SM process qq — > g — > ti, 
the ti pair can also be produced via a t-channel process with a W mediator. Apart from a 
common factor —ig^gp^E^/ (t — rri^,), the helicity amplitude M^,(A g , Xq, X t , Af) is 

MU+ ) = - [2 + r 2 w ] T^T^sine 

= [2(l-0)+r&(l + /3)](l-cos0) 

M\y,(+-+-) = [2(l + 0)+r&(l-0)](l + cos0) 

M t w ,(+-++) = [2 + r 2 w ] Vl-/3 2 sin^, (25) 

where rw = m t /mw'- 

In the region /3 ~ 1, the nonzero helicity amplitudes are 

M\y,(+--+) ~ 2r^(l-cos^), 

M t w ,(+ - +-) ~ 4(1 + cos^) . (26) 
12 



In order to produce top quarks in the forward region, one needs 2r^ < 4, which is always 
true for the region of W masses (heavier than the top quark) considered in this paper. At 
the Tevatron the ft distribution of the top quark in tt production peaks around 0.6, and 
therefore most of the top quarks are not significantly boosted. We can also easily see that 
more right-handed top quarks are produced compared to left-handed ones in the W model, 
p tR > pt L . Since the t-channel propagator contributes a minus sign, the total forward- 
backward asymmetry results from a competition between the square of the purely NP term 
and the interference term of NP with the SM. The former is proportional to g\ and the 
latter to g\. We plot the correlation between A FB and A e FB for the W model in Fig. IU(b). 
The strong correlation is fit well by 

A FB ~ 0.75 x A FB - 2.1% . (27) 

Moreover, for a relatively light W 600) GeV, both A FB and A FB can be consistent with 
the DO data within 1 a. 

For the G' and W models, Fig. H^c) shows the statistics for the ratio of the predicted A FB 
to A FB , based on the scattered points in Fig. Ufa) and H£b). The total number of events is 
normalized to 1. The axigluon model peaks near 50% and W model near 62%. The ratio in 
the SM is close to 40%. The W model generates a larger A FB than the axigluon G' model 
because it produces more right-handed top quarks. The comparison to the DO point shown 
in Figs. ID^a and b) indicates that top quark events with a large proportion of right-handed 
top quarks are favored. Constraints on flavor-changing currents in the W model allow only 
right-handed couplings to the top quark, consistent with the DO A l FB results. There is no 
direct evidence of the handedness of the coupling in the massive gluon models. The DO 
result could be interpreted as an indirect clue for the chiral couplings of the massive gluon. 
Improved statistics would help, as well as a measurement of A FB by the CDF collaboration. 

V. CONCLUSION 

The deviation of the top quark forward-backward asymmetry A FB from its SM prediction 
may indicate the presence of new physics. Based simply on the large value of A FB , the 
charged lepton forward-backward asymmetry A l FB should also be expected to be larger 
than the SM expectation. Indeed, the DO collaboration reports A FB = 15.2%, about 3a 
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away above the SM value. In this paper, we study the kinematic and dynamic aspects 
of the relationship between the asymmetries A FB and A l FB based on the spin correlation 
between charged leptons and the top quark with different polarization states. Owing to 
the spin correlation in top quark decay, A FB and A FB are strongly positively correlated for 
right-handed top quarks. However, for left-handed top quarks, the nature of the correlation 
depends on how boosted the top quark is. For large enough top quark energy, left-handed top 
quarks will also generate a large charged-lepton asymmetry, similar to that for right-handed 
quarks. However, if the top quark is not boosted (E t < 200 GeV), A FB from left-handed top 
quarks will be less than A FB /2 for a positive A FB . Since most of the ti events are produced 
in the threshold region, one may use the large positive values of A FB and A FB measured at 
DO to conclude that production of left-handed top quarks is disfavored. Confirmation of the 
DO result and greater statistics are desirable. There is great value in making measurements 
of both A FB and A l FB because their correlation can be related through top quark polarization 
to the underlying dynamics of top quark production. 

We focus on two benchmark NP models, an axigluon (G") model which produces un- 
polarized top quarks, i.e. an equal number of right-handed and left-handed top quarks, 
and a flavor-changing W' model which produced dominantly right-handed top quarks. To 
determine free parameters, we require that these new physics models fit A FB as well as 
the tt total cross section at the Tevatron at la level of accuracy. As we show, there is a 
strong correlation between A FB and A FB in both models. The best fit to the relationship 
is A FB ~ 0.47 x A FB + 0.25% and A FB ~ 0.75 x A FB - 2.1%, for the axigluon model and 
the W model, respectively, both within 2a of the DO result. To generate A FB satisfying the 
data to better than la accuracy, a heavy G' (heavier than about 1 TeV) is preferred, and a 
light W (lighter than 600 GeV) is favored. 

We do not address the LHC case in this paper but may do so at a later time. Owing 
to the lack of definition of a forward direction in a pp collision, it is less straightforward to 
measure the two observables we discuss here. 

Acknowledgments 

The work of E.L.B., C.R.C. and H.Z. is supported in part by the U.S. DOE under 
Grants No. DE-AC02-06CH11357. H.Z. is also supported by DOE under the Grant No. 



14 



DE-FG02-94ER40840. The work of J.H.Y. is supported in part by the U.S. National Science 
Foundation under Grand No. PHY-0855561. 



[1] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. D 83, 112003 (2011) 
[2] V. M. Abazov et al. [DO Collaboration], la£XiyJ407^995] [hep-ex]. 

[3] J. H. Kiihn and G. Rodrigo, Phys. Rev. Lett. 81, 49 (1998); J. H. Kuhn and G. Rodrigo, 

Phys. Rev. D 59, 054017 (1999). 
[4] D. Krohn, T. Liu, J. Shelton and L. T. Wang, arXiv:1105.3743l [hep-ph]. 
[5] A. Falkowski, G. Perez and M. Schmaltz. larXiv: 11 10.37961 [hep-ph]. 
[6] S. Jung, H. Murayama, A. Pierce, and J. D. Wells, Phys. Rev. D81, 015004 (2010). 
[7] K. Cheung, W. Y. Keung and T. C. Yuan, Phys. Lett. B 682, 287 (2009); 
[8] P. Ferrario and G. Rodrigo, Phys. Rev. D 80 (2009) 051701; 
[9] P. H. Frampton, J. Shu and K. Wang, Phys. Lett. B 683, 294 (2010); 
[10] Q.-H. Cao, D. McKeen, J. L. Rosner, G. Shaughnessy, and C. E. M. Wagner, Phys. Rev. D81, 

114004 (2010). 

[11] B. Xiao, Y. -k. Wang, S. -h. Zhu, Phys. Rev. D82, 034026 (2010). [arXiv:1006.2"5T0l [hep-ph]]. 
[12] K. Yan, J. Wang, D. Y. Shao and C. S. Li. larXiv:1110.6684l [hep-ph], 

[13] E. L. Berger, Q. H. Cao, C. R. Chen, C. S. Li and H. Zhang, Phys. Rev. Lett. 106 (2011) 
201801 

[14] S. Chatrchyan et al. [CMS Collaboration], JHEP 1108 (2011) 005 
[15] G. Mahlon and S. J. Parke, Phys. Rev. D 53, 4886 (1996) 

[16] E. L. Berger, Q. H. Cao, C. R. Chen and H. Zhang, Phys. Rev. D 83, 114026 (2011) 



15 



